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Quantum random walk of two photons in separable and entangled state 
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We discuss quantum random walk of two photons using linear optical elements. We analyze the 
quantum random walk using photons in a variety of quantum states including entangled states. We 
find that for photons initially in separable Fock states, the final state is entangled. For polarization 
entangled photons produced by type II downconverter, we calculate the joint probability of detecting 
two photons at a given site. We show the remarkable dependence of the two photon detection 
probability on the quantum nature of the state. In order to understand the quantum random walk, 
we present exact analytical results for small number of steps like five. We present in details numerical 
results for a number of cases and supplement the numerical results with asymptotic analytical results. 

O ' PACS numbers: 03.67.Lx, 42. 25. Hz 
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^ ; I. INTRODUCTION 

^^ . A new_paradigm in the study of random walks has recently emerged [ll, [2, S 0, B H 0, H M, llO, [ill, E, [iJ, H, 

llSl . [la . [l7| . The random walker has been assigned an additional quantum degree of freedom which could be spin 

degree of freedom [1| . Thus walker goes left or right depending on the spin degree of freedom. The probability of 

^f^ • finding the walker at a given site now depends on the spin state of the walker. All this has now been well studied 

^ ! [l| and numerical simulations have been done to find the site distribution after walker has taken large number of 

OO . steps. Using linear optical elements, Do et al [g| have realized the quantum random walk (QRW). However, in their 

^^ ■ experiments they used a weak coherent field rather than a field with strong quantum character. This is fine if we 

recall the results of Knight et al [ll| who showed that QRW of a single walker can be implemented by using classical 

fields. A similar arrangement is discussed earlier by other authors as well [3, [l3|- Jeong et al [10| analyzed the cases 

^^ , of a walker in a coherent state and in a single photon state and concluded that the final probability distribution was 

f^ ■ identical in the two cases; though different from that of the classical random walk (i.e. a walker without the additional 

quantum degree of freedom). This is explained by Jeong et al \10\ in terms of the P-representation of the state of 

photons. Thus an important question is-what would be a strict QRW which can not be produced by using classical 

fields. To understand this aspect, we study QRW by two photons in a variety of quantum states including entangled 

"S ,. states. We find that even if initially the two photons are in separable Fock states, the final state is entangled. This is 

C^ ' quite an interesting quantum property and has no classical counter part. Omar et al [l6| have studied QRW of two 

^ nonidentical walkers with entangled initial state and have shown that the final state is entangled depending on the 

O"^' initial entanglement. Our model of QRW is different from that of Ref [1^ as our system can produce entanglement 

►>. I even if initially there is none. This is clearly borne out by our result in Sec. IV A. Further, we calculate photon-photon 

• ""j ■ correlations which in the past have been very successfully used to reveal the quantum character of the fields [la, [l^ . 

KN ' We show the remarkable dependence of the two photon detection probability on the quantum nature of the state of 

^ . photons. We present explicit results for all four Bell states of the incoming photons. Our work thus brings out the 

■ - - ' role of coherences and entanglement in QRW. 

The organization of the paper is as follows. In Sec. II, we present the system of linear optical elements used to 
realize QRW. In Sec. Ill, we discuss QRW of a single photon when the photon enters in the arrangement through 
the two input ports. The state of the incoming photon is a pure state. We present an approximate analysis for the 
probability of finding the photon at a site after a large number of steps |2d . [2l| and compare with the exact numerical 
results. In Sees. IV and V, we consider QRW of two photons. We derive analytical results for the photon-photon 
entanglement corresponding to the case when the walkers take only a small number of steps. The analytical results 
clearly demonstrate the dependence of two photon detection probability on entanglement between the walkers. We 
evolve a numerical strategy as well as an approximate analysis to obtain results for final state of the photons after a 
large number of steps. In Sec. IV, we consider the case when the initial state of the photons is separable however the 
final state of the photons is entangled. Here the entanglement is generated due to the passage of photons through 
linear optical elements. We also discuss a case when two photons in Fock state are replaced by two photons in coherent 
states. We point out that QRW in our scheme with photons either in separable Fock states or in an entangled state 
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FIG. 1: The schematic arrangements for reahzing QRW of two entangled walkers. In the inset effect of polarization beam 
splitter is shown explicitly. The x-polarized photon is transmitted in the same direction but the y-polarized photon is reflected 
and changed the direction of propagation. 

can not be reproduced by using coherent states. In Sec.V, we discuss QRW of two photons when initial state is an 
entangled state and discuss the dependence of two photon detection probability on entanglement initially present in 
the state of the walkers. We present our conclusions in Sec. VI. 



II. DESCRIPTION OF THE ARRANGEMENT 



In Fig. [Tl we show a schematic arrangement for realization of QRW on a line using linear optical elements like 
polarizing beam splitters (PBS) and half wave plates (HWP). The arrangement looks like a large interferometer and 
has been discussed in an experimental realization of quantum quincunx [8[ . The photons enter in the arrangement 
through the two input ports. The state of the photons anywhere along the arrangement can be completely defined in 
terms of its direction of propagation and state of polarization. A half wave plate performs Hadamard operation for a 
single photon in polarization state basis [23 |. 
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where \x) and \y) are two orthogonal linear polarization states of the photons. The polarizing beam splitter (PBS) 
produces polarization dependent spatial displacement in the position of the photons. Clearly the displacement will 
be either in horizontal or in vertical direction depending on the direction of propagation and the state of polarization 
of the incoming photons. The operation of PBS can be written as 
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where h{v) shows the propagation of the photons along the horizontal (vertical) direction. After passing through PBS 
photons are displaced by one step in the horizontal direction when they come in the state \h, x) or jw, y) and displaced 
by one step in vertical direction when they come in the state \h, y) or |w, x). Thus, in our scheme, states |/i, a;), \h, y), 
\v,x), and \v,y) are corresponding to the states of a four sided coin and the position of the photon along the optical 
arrangement performs QRW. 

The photons enter through the input ports and pass through a half wave plate followed by a PBS. The half wave 
plate performs Hadamard operation in the polarization space. Then PBS produces shift in the position of the photons 
depending on the polarization and the direction of the propagation of the incoming photons. Thus each layer of half 



wave plate followed by the PBS generates one step of QRW in the spatial state of the photons. After particular 
number of such steps the position of the photons is detected by an array of the detectors. The detection of the 
photons can be made polarization sensitive or insensitive depending on the particular interest of measurement of the 
state. 

In Fig. [U we show the arrangement which has five layers, each layer is made of half wave plate followed by a PBS. 
This arrangement can produce QRW by five steps. Needless to say that for large number of steps a large number of 
optical elements will be required and such a big arrangement may have some difficulties in alignment as well as the 
visibility of the signal. 

III. QRW OF A SINGLE PHOTON 

In this section, we consider the case when a single photon enters through the input ports. Thus we have a single 
walker and the QRW is controlled by a four sided coin. The four sides of the coin are formed by the four degrees of 
freedom of the photon, two directions of propagation as the photon can move in horizontal or in the vertical direction, 
and two directions of polarization. The state of the photon at any time during the QRW can be defined as |n, g, j, fc), 
where n represents the number of steps, q represents the spatial position of the photon with the convention that the 
horizontal displacement is in -\-ve direction and the vertical displacement is in —ve direction, and j, k represent the 
direction of propagation and the direction of polarization respectively. 

The dynamics of the photon can be described by expressing the state of the photon after n + 1 steps in terms of 
the states after n steps in the following way. 
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First the half wave plate generates transformation in the polarization states of the photon which is followed by 
the spatial transformation generated by the PBS. We can reexpress Eqs. (131)-® in terms of operators C and S, 
equivalently, as 
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where C and S are 



C= ^ i\h,x){h,x\ + \h,x){h,y\ + \h,y)(v,x\ - \h,y){v,y\ 
+ \v,x){v,x\ + \v,x){v,y\ + \v,y){h,x\ - \v,y){h,y\) 



(8) 



S=\h,x){h,x\ ^ \q+l){q\ + \h,y){h,y\ O \q+l){q\ 
+ \v,x){v,x\ <»\q- l){q\ + \v,y){v,y\ (g) |g - l){q\ 
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FIG. 2: (Color online) The probability distribution P(q) for detecting the walker at position q, for the initial coin state (a) 
-y={\h,x) + \v,y)) and {h)^{\h,x exact (approxiinated) values. 
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FIG. 3: (Color online) The probability distribution P{q) for detecting the walker at position q, for the initial coin state (a) 
-y={\h, y) -f \v, x)) and (b)-i=(|/i, y) — \v, x)). The black solid (red dashed) line is corresponding to exact (approximated) values. 



Thus, the unitary coin operator C and the conditioned shift operator S for the QRW of the photon are given by 
Eqs. ([5]) and (O respectively. The transformation Ug = S (E) C acting on the initial state of the photon is equivalent 
to the one step of QRW. The iterative application of the transformation Us for n times, C/", gives the QRW of n 
steps. The initial state of the coin can be selected as one of the states, \h,x), \h,y), \v,x), \v,y), and their possible 
superpositions. For a classical analog to this QRW we look at the shift operator S defined by Eq.([S]). Clearly, this 
walk is equivalent to the random walk of a single walker with four dimensional coin where walker moves one step 
forward for two possible outcomes of the coin and moves one step backward for the other two. The probability of 
reaching the walker at position q will be a binomial distribution 
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Thus classically we can not differentiate between the walk of single photon in this case and the walk of a single walker 
with two sided coin. In the following we show how QRW of single photon in this case leads to different results than 
QRW with a two sided coin. 

Here we present an approximate analysis of QRW for large number of steps. We follow a similar method as discussed 
by Nayak et al [20] and Brun et al [2l| using discrete spatial Fourier transform and taking asymptotic limit of the 
Fourier integrals for large number of steps. We relegate the details to the appendix. For initial coin state \h,x), we 



get the final state of the photon |^/ia;(n)) after n steps as follows. 

n 
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q— — n 

ItphxiiT^q)) == fhx\n,q,h,x) + fhy\n,q,h,y) + fyx\n,q,v,x) + fvy\n,q,v,y), 
where the coefficients /^^ are given by 
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where (f){k) = —{tOk + ka), a = q/n, ujk € [7r/4, 37r/4] and defined as coswfc = cosfc/-\/2, fco = sin~"^(— a/Vl — a^) 
and prime denotes the derivative with respect to k. This approximation is vahd in the interval a G [— l/v2, 1/ v2], 
outside of this interval f^y can be taken zero. In the expressions of /^j^, the integral inside the bracket is independent 
of n, and is responsible for the constant spikes in the probability of detecting photon near initial position g = 0. We 
evaluate this integral numerically. The second term inside the bracket is sum of cosines and completely characterizes 
QRW. Similarly, for initial coin state |/i, y), the state of the photon after n steps is 
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From the symmetry of the arrangement (see Fig. [T]), the state of the photon after n steps for initial coin states \v,x) 
and \v,y) can be written by interchanging the direction of propagation h and v and replacing q by —q in \^hx{n)) 
and \'^hy{n)) respectively. 
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For initial coin state as an arbitrary superposition a\h,x) + f3\v,y), the state of the photon after n steps is given by 

|*H) = a|*ft,(n)) + Pl^.yin)). (28) 

In Figs. [2] and [HI we have shown the probabihty distribution for QRW of the single photon with initial coin states 
-y={\h,x) ± \v,y)) and -y={\h,y) ± \v,x)) after 50 steps. We plot the results using above approximate analysis as well 
as exact numerical simulations. Clearly for large number of steps, say for n = 50, there are very small differences 
between the approximate analysis and the exact simulations. Further for larger values of n these differences will be 
negligible. It should be noted that the probability distributions for detecting the photon at position q in Figs. [2] 
and [3] are very much different than the distributions for QRW with a two sided coin Q . In all these cases the initial 
state is most probable state and the distribution is sharply peaked at q — 0. For initial coin states -y={\h,x) + \v,y)) 

and -j=(^\h,y) + \v,x)) the distributions are symmetric but the side peaks are very small and most of the time walker 
remains at its initial position very precisely. In the case of QRW with initial coin states -j={\h,x) — \v,y)) and 
-j^{\h,y) — \v,x)), in addition to the narrow central peak, distributions have a peak along one side of the position 
axis. Further for the state -j^{\h,x) — \v,y)) the additional peak in the distribution is along the positive side of the 
axis at g = n/^/2, while for state -n={\h, y) — \v, x)) the additional peak is along the negative side a,i q — —n/\/2. We 
emphasize that the quantum random walk of a single photon depends very much on the initial state, see for example 
the distinction between the Figs. [DJa) and^b). 

IV. QRW OF TWO PHOTONS WITH SEPARABLE INITIAL STATE 

In recent papers [ll|, it has been shown that QRW is an interference phenomenon and does not essentially depend 
on the quantum nature of the state of the walker. As a result various classical sources like low intensity lasers [^, [iJl 
and coherent state of radiation fields 0, [lj| are used to realize QRW. In order to explore further the quantum nature 
of random walk, we consider the case when two photons start QRW from a separable initial state 

|*) = |*i)®|*2), (29) 

where |$) is state of two photons and |^i) and |^2) are states of single photons. The state of the photons after a 
certain number of steps is, in general, not a separable state as quantum entanglement is produced by linear optical 
elements. This is reminiscent of the well known property [l8[ of a 50-50 beam splitter where two incoming photons 
in the separable state |1, 1) go over to an entangled state of the form (|2, 0) + |0, 2))/V2. We first consider the case 
of input states which are single photon states. We would also consider the case when the input states are replaced by 
coherent states 

A. QRW of two photons with initially in separable Fock states 

In our scheme, two photons act as two walkers. They enter in the arrangement through the two input ports. Initially 
one photon propagates in horizontal direction and the other in vertical direction. We consider the initial state of the 
photons as one of the four separable states |0,0, h,x) ® |0, 0,w,x), |0,0, h,y) ® |0,0, w,y), |0,0, h,x) ® |0,0, w,y), and 
|0, 0, h, y) (g) |0, 0, w, x). We can write these states in terms of initial field operators as follow. 

|0, 0, h, x) ® |0, 0, V, x) EE ai*^°^^air'^ |0), (30) 

|0, 0, h, x) ® |0, 0, V, y) EE afj'^^a^^f^ |0), (31) 

|0, 0, h, y) ® |0, 0, V, y) ^ ai^/^^ai^^t |o), (32) 

|0,0,/i,2/) ® W\v,x) ^ ai^^^air^^io), (33) 

where a,^''' is creation operator for a photon, a,'^'* |0) = \n,q^i,k) ^ and |0) is vacuum. Here we present the 
analytical calculations for QRW of few steps, say five steps. For linear optical elements, it is sometimes instructive 
and transparent to work with the transformation of operators. This is particularly so if the states with more than 
one photons are involved. Thus for calculating final state of the walkers after five steps, we express the initial field 



operators, aL' , in terms of the final field operators, a.ji^'^ , after five steps. 
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Using transformations ([M]) to ([57]) and Eqs. ((5D|) to ((55|) . the state of the quantum walkers, corresponding to a 
particular initial state, after five steps can be calculated. Clearly, the final state of the photons is an entangled state 
and can not be expressed as a product of the states of two single photons. 

I*) 7^1*1)® 1*2), (38) 

where |^) is final state of two photons and |^i) and [^2) are final states of single photons. It should be borne in 
mind that Eqs. ([34]) to ((37)) should be supplemented by free field operators at the open ports. These are important for 
the operator algebra. However these do not contribute to the results below and hence for brevity we have not written 
these explicitly in Eqs.dMl) - dSB- 

Finally, the results for calculated probability P{qi, 92) for detecting the walkers at positions qi and 92 simultaneously, 
after five steps, are shown in the Table. I. Note that the diagonal elements P{q, q) give the probability of finding two 
photons at the site q. It is clear from the Table. I that the probability distributions for all considered initial states are 
different to each other. At the bottom of the table we present the probability P{qi) of detecting at least one photon at 
the position qi, where P{qi) — J2q2 PiQi^l^)- Further using the values of P(gi, 92) and P{qi), P{q2) for the positions 
qi and 52 we can also calculate the correlation 

(Tq,c„^P{qi,q2)-P{qi)P{q2)- (39) 

We found that the correlation Uq^q^ is non zero almost everywhere for all values of P{qi,q2), which shows that 
though initially walkers were in a separable state, but after few steps their state is entangled. We emphasize that the 
correlation ([39]) is due to quantum as the state (p8)) does not factorize. These correlations are due to the quantum 
nature of the initial state and arise when the photons pass through the linear optical elements. The origin of such 
correlations has been observed by Mandel and coworkers '13| in their pioneer work on beam splitters. It should be 
noted here that if we replace the photons with two coherent states the output state will be a factorized state and the 
probability distribution will not exhibit such correlations. Thus the QRW of two photons is completely dependent of 
quantum nature of the state of the photons and no coherent state can reproduce such QRW. It should be noted that 
the normalization condition for P{qi,q2) is given by ^ ^ P{qi,q2) — 1. Further ^ P{qi) is not equal to 1 and 
the normalization condition for P{qi) will be ^ P{qi) + P{qii 9i) — 2. To see it more clearly consider the following 
state of finding two photons at sites 1 and 2, 

IV') = ^(|1,1) + |1,2) + |2,2)) (40) 

For this state the probabilities of detecting at least one photon at site 1 and 2 are P(l) — -P(2) — 2/3 and the 
probabilities of detecting both photons at same site are P(l, 1) = ^(2, 2) = 1/3 which satisfy the above normalization 
condition. 



TABLE I: The calculated probabilities of detection for the photons after five steps, for initial state (a) \0,0,h,x) ® \0,0,v,x), 
(b) jo, 0, h, x) ® |0, 0, V, y), (c) |0, 0, h, y) ® |0, 0, v, y), and (d) |0, 0, h, y) ® |0, 0, v, x). All values shown in the table are 128 times 
of the actual values. 
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FIG. 4: (Color online) The probability P{qi,q2) of detecting the walkers at positions 91 and 92 after number of steps n = 25. 
The walkers come in the separable state (a) \0,0,h,x) (g) |0, 0,w,a:), (b) \0,0,h,x) <Si \0,0,v,y), (c) \0,0, h,y) (g) \0,0,v,y), and 
(d) \0,0,h,y)<S,\0,0,v,x). 



After providing an approach to the analytical calculations for few steps, we present numerical simulations for larger 
number of steps. The operator U = (S* (8) C) (8) (5 (8) C) acting on the initial state of the photons generates QRW of 
one step. Thus the transformation U" will give the final state of the random walkers after n steps. Here unitary coin 
operator C and the conditioned shift operator S for QRW are given by the Eqs. ^ and ([9]) of the previous section. 

In Fig. Uwe have plotted the probability distributions for detecting photons at positions qi and 52 simultaneously 
after 25 steps of QRW. In Fig. |3]the photons are in separable initial states. We notice that each probability distribution 
is very much different from the other. A very common feature in all plots is sharp central peak which shows that in 
all cases the initial state is the most probable state. In fact this is the property of classical random walk, which has 
Gaussian probability distribution, but in Fig. 2] the central peak is much narrower than a Gaussian distribution and 
the large spread of the distribution has strickenly different behavior. 

For approximate analytical calculations for large number of steps we use the results (|12p , ((T9| , ((26)) and ([27)1 in the 
following way. For example, for initial state |0, 0, h, x) (g) |0, 0, v, y), the state of the photons after n steps will be 



l*xs,(ri)) 



1 



V2. 



y^ [\->Phxin, qi)) (g) \ipvvin, 92)) + \iphxin, 92)) ® {-ipvyin, qi))] 



(41) 



Here I'ixyin)) has exchange degeneracy as both photons are undistinguished after they have passed through the 
optical arrangement. Similarly, we can calculate the final state of the photons for other separable initial states shown 
in Figdl We have checked that the joint probability distributions P{qi, (72) for large number of steps calculated using 
above approximate analysis match with the exact simulations. Thus our approximate analysis can be used to calculate 
the correlations between the walkers after large number of steps very well. 

Here we have explicitly shown that in our case QRW of two photons is highly entangled, even though the photons 
start in separable states. The entanglement is developed in the course of time when the photons pass through the 
optical arrangement. In the next subsection we show that the QRW of two photons in our case can not be produced 
by using coherent states. 



TABLE II: The normalized probability P(q) of detecting single photon at particular position q, after QRW of 5 steps with 
initial state \a)f^^ \P)vy , for small amplitudes a and /3. 



a 


/3 


Probability P(g) 
q = —5q = —3q = —lq = l q = 3 q = 5 


f 
f 
f 



/ 

-/ 


/ 


1/16 1/8 3/8 3/8 1/16 

1/16 1/8 3/8 3/8 1/16 
1/32 3/32 3/8 1/8 11/32 1/32 
1/32 3/32 1/8 5/8 3/32 1/32 



B. QRW of two photons initially in separable coherent states 



Here we discuss the case of QRW when two photons in our scheme are replaced with two weak coherent states. Let 
us consider that the initial state is 



1^(0)) 



a 



•r'l 



\(0,0) 

/ vy ' 



(42) 



where indices to the coherent states let) and 
field operators can be expressed as |24| 



have their earlier assigned meanings. The initial state in terms of 



IV'(O)) = expiaa^^^ ~ a^a^"/)) exp(/3ai°/)t - /?*a(°/^)|0). 
Now using transformations ([34|), (I37p and Ea. (|43| the final state, after five steps of QRW, is given by 



(43) 



|V(5)) 



I " + P x(5.-5) I Q^ + /^ N (5,-3) I a - P ,(5.- 

'4V2^'"" UV2'''^ '4^/2^"" ^yz 

^| Q' + /^ ^(5,-i)| 3a + /3 ^(5,-i), a + /3 ^(5 _i). a~ P 
' 4^2 ^''" ' 4^2 '''y ' 4^/2 ^"" ' 4^/2 



(5,-3) 
vy 



a — P , 

-J-ff/v. l^-^/hv l^^/v. I^yl-. 

-VZ 4VZ ^ 4V2^"" '4V2^"^ 

| Q + 3/3 (5,i), Q-3/? ^(5,i)| Q-/3 ^(5,i), -a-/? ^(5,i) 
'' 4^2 ^'" ' AV2 ^"'^ ' 4V2 ^"" ' AV2 ^"^ 

-3a + /3 ^(5.3)| "« + /3 ^(5.3)| a-/3 ^(5 3), Q-/3 (5,5) 

' 4^/2 '''^ ' 4^/2 ''^y ' 4V2 '^^ ' 4^2 ^''" ' 



(44) 



where superscript (5, q) on each state in (|H)l denotes the position of the photons at q after 5 steps. Clearly the final 
state (|M1) is a product state and the value of correlation aq^q^ defined by ([55)1 will be zero for such state. 

Note that the probability of finding one photon in a coherent state ja) is jape" I"' . For small enough \a\, it reduces 
to jap which is the mean number of photons in a coherent state. Thus, the normalized probability P{q) of detecting 
a photon at q in this case is equal to the average number of photons detected at site q divided by the average number 
of incident photons |ap + |/3p. 

In Table. II, we show the normalized probabilities of detecting photon at site q for small values of a and /?. By 
comparison of the cases corresponding to a = /, /3 = / and a — f, f3 — 0, where |/| <C 1, we see how the interference 
effects change these probabilities. Similarly for a — f, /3 — f and a = 0, /3 = f, the probabilities depend on the 
interference among various paths. Further these probabilities for one photon detection are different from those listed 
in the Table. 1(b), for the case of two photons in separable Fock states \0,Q,h,x) (g) \0,0,v,y). It is noted that the 
normalized probability distribution of detecting photons at particular position q for such initial states is identical to 
the case of single walker discussed in Sec. Ill, not to the case of two walkers in Sec. IV A. For example, for a = f and 
P = f the distribution is identical to Figl^l^a) and ior a — f, P = —f we get distribution identical to Figj^JJb). This 
can be explained by examining the projection of state (|42p in one photon space, |-0o) — > Oi\h, x) ® |0) + /3|0) |w, y) as 
the contributions from spaces containing more than one photon are of higher order in a and p. 



V. QRW OF TWO PHOTONS IN ENTANGLED STATE 



Here we consider the case when photons start QRW in an entangled state. Thus one photon enters from each 
input port and the state of the photons is maximally entangled in polarization basis, as those produced by a type-II 
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TABLE III: The calculateid probabilities of (detection for the photons after five steps, for initial state (a) 

1/^2 (|0, 0, h, x)(g)\0, 0, V, y) + |0, 0, h, y) ®\Q, 0, «, x)), (b) 1/a/2 (|0, 0, h, x) ® |0, 0, v, y) - |0, 0, h, y)®\0, 0, v, x)), (c) 

l/^/2(|0,0,/l,a;)(^|0,0,^;,a;) + |0,0,/l,y)(^|0,0,w,y)), and (d) 1/^2 (|0,0, /i,a;) ® |0,0, w,2:) - |0, 0, /i, y) ® |0, 0, v,?/)). All 
values shown in the table are 128 times of the actual values. 
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FIG. 5: (Color online) The probability P{qi,q2) of detecting the walkers at positions 91 and 92 after number of steps n = 25. 
The walkers come in the entangled state (a) I'tp'^), (b) Itp^), (c) \(l)'^), and (d) \(t>^). 



downconverter. We consider that the state of the photons at input ports is one of the four Bell's states (2E 

IV-*) = -^(|0, 0, h, x) ® |0, 0, V, y) ± |0, 0, h, y) ® |0, 0, v, x)), 
|(^±) = -^(|0, 0, h, x) ® |0, 0, V, x) ± |0, 0, h, y) ® |0, 0, v, y)). 



The states (|45| and (|46| can be written, in terms of the field operators, as 



(45) 
(46) 

(47) 
(48) 



Following a similar procedure as discussed in Sec.IV, we can calculate the final state of photons after five steps using 
transformations ([M)) to ([57]) and Eqs. (|47p to (|^5)) . In this case also we find that the final state of the photons is 
an entangled state. Thus in this case both initial and final states of the walkers are entangled state and no classical 
analog for such states is possible. 

In Table. Ill, we show the results for detecting the photons simultaneously when they start QRW from one of 
the Bell's states (|47l) and P5|) . We notice that the differences between the distributions for symmetric states and 
antisymmetric states are much larger than the differences between the distributions for symmetric-symmetric or 
antisymmetric- antisymmetric states. In terms of correlation (|39p . both the initial and the final states are highly 
correlated. 

For larger number of steps, we do numerical simulations following the method discussed in Sec.IV A. In Fig. [S] 
the initial states of the photons are maximally entangled states (P5|) and (|^5|) . For different entangled states we get 
different probability distributions. The differences between the distributions for symmetric states and antisymmetric 
states are much larger than the differences between the distributions for symmetric-symmetric or antisymmetric- 
antisymmetric states. Further, we have seen by direct computation that the probability distributions for entangled 
states can not be reproduced exactly by the walkers having initial states as an incoherent mixture of separable states. 
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FIG. 6: The probability distribution P{qi) for detecting at least one walker at position qi after 100 steps. For initial states 
\tp ) (dotted line) and \4>^) (solid line), the distribution remains almost same. 

In Fig. [6] we have shown the probability distribution P{qi) for detecting at least one walker at the position qi. 
In this case the distribution is symmetric for both negative and positive values of qi and has a sharp maxima at 
initial position qi — 0. The probability distribution for walkers with maximally entangled initial states (j45p and (j46p 
remains almost invariant except at initial position. At initial position the probability of detecting one walker is larger 
for states (|45|) than the case of initial state ([46]) . In the case of two-photon QRW with separable Fock states discussed 
in Sec. IV A, the probability P{qi) is different for different initial state. For initial states \0,0,h,x) (E) \0,0,v,x) and 
\0,0, h,y) (8) \0,0,v,y) the distribution P{qi) is symmetric on both sides of qi-axis and side peaks appears in both 
positive and negative directions, while in the case of |0, 0, h, x) ® |0, 0, u, y) and |0, 0, h, y) ® |0, 0, w, x) the probability 
distributions are antisymmetric with one side peak in positive and negative direction respectively. 

For large number of steps, we can again use our approximate results P^ . p^ . (^5)1 and ([?71) in the following way. 
For initial state as one of the Bell's state |^^) the final state of the photons after n steps will be 



|vl/±(n)) 



E 

91,12=-" 



[\%ljhx{n,qi)) ® \^yy{n,q2)) + \iphxin,q2)) » |'0i,j^(n, gi)) 

±\i>hy{ri,qi)) (g) \'ipvx{n,q2)) ± \i>hy{n,q2)) (gi |?/'t,x('^, <?i))] • 
Similarly, for initial state as Bell's states \<j>^), the state of the photons will be 



(49) 



1 

91,92=-" 



hx{n,qi)) «) \ipvx{n,q2)) + \'iphx{n,q2)) 'S) \'ipvxin,qi)) 
±\'^hy{n,qi)) (g) \ipvv{n, q2)) ± \iphy{n, q2)) <E) \i(;vy{n,qi))] . 



(50) 



The results ([^^ and ([50)1 match very well with exact simulations for large number of steps. Here it should be noted 
that QRW in our case is different than the case of two entangled walkers discussed in Ref. [16|]. In Ref. T6], the 
correlation between the walkers is completely because of their correlated initial state only. The walkers are nonidentical 
and evolve completely independent to each other. Thus if the initial state will be a separable state there will be no 
correlation in QRW and can be produced using walkers in coherent states. It should be borne in mind that we can 
not construct entangled state of one x-polarized photon and the other y-polarized photon unless some other degree 
of freedom like direction of propagation is introduced. 



VI. CONCLUSIONS 



In conclusion, we have discussed how the quantum features of QRW can be uncovered by studying QRW by 
photons in a number of quantum states and by detecting the coincidence correlations between two photons. We use 
two photons, each of which can exist in two polarization states. Further, the photons can travel in either vertical or 
horizontal direction. The walk is realized by optical elements consisting of polarization beam splitters and half wave 
plates. We first consider the QRW by a single photon and show how the result of QRW depends on initial state of 
the photon. A comparison of the Figsj^fa) and[2I^b) shows that the probabilities depend on the relative phases in the 
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initial superposition state of the photon. We next consider QRW by two walkers in separable quantum states like 
Fock states. We show that our optical arrangement entangles two photons even though initially they are in separable 
states. The joint probability of finding the two photons at two different sites now depends on the initial separable 
quantum state. Further the probability for single photon detection is different from those calculated in Sec. III. We 
also examine the QRW of two photons in coherent states. In this case the state of the two photons remains separable. 
The single photon detection probabilities depend on the initial amplitudes of the two coherent states and thus the 
interferences are prominent. Finally we consider QRW of two photons in all four Bell states. The resulting joint 
probability distributions are quite different from those calculated for separable states. All the cases discussed reveal 
very interesting quantum character. Clearly our analysis is applicable to other systems as well, for example electrons 
with appropriate arrangement of Stern-Gerlach fields. 
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APPENDIX A: FOURIER ANALYSIS OF QRW OF A SINGLE PHOTON 

We can write Eq. ([7|) for the wave function of the single photon in the form 

V'(n + 1, g) = M+ij{n, q-l)+ M-^{n, q + 1) 
where 



"^ - 71 



V2 



110 0" 
1-1 

0. 
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(Al) 



(A2) 



(A3) 



and the state of the photon after n steps, \'4>{'n, q)) = Chx\n, q, h, x) + Chy\n, q, h, y) + c^xW, q, v, x) + Cyy\n, q, v, y), is 
expressed in the matrix form 



V'(n,g) 



Chx{n,q) 
Chy{n,q) 
Cvx{n,q) 
Cyy{n,q) 



(A4) 



Now we solve Eg. ljAip using spatial discrete Fourier transform. The spatial discrete Fourier transform for k e [— tt, tt] 
is defined by 



^(fc) = ^V(x)e''^ 



and ipix) 



27r 



tP{k)e-"'='dk 



(A5) 
(A6) 



We write Eq.jATJin the Fourier domain using transform (jASP as following, 

^(n+l,fc) = ^M+i!{n,x-l)e'^^ +M^ij{n,x + l)e^^'' 

X 

= {e^M+ + e-^^MJ)i){n, k) 
= Mki>{n,k) 



(A7) 
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where 



Mfe = 



V2 



g,:fc 


^ik 





1 








^ik 


_gifc 








^-^k 


g—ik 


^-ik 


_g-»fc 









(A8) 



It is clear from Eq. (jA7p that in the Fourier domain the n step QRW of the photon starting from the state -(/"(O, k) is 
given by 



Vi(n,fc) = (A4)"V;(0,fc) 

= ^(A™)"|*™)($™|^(0,A:)) 



(A9) 



where A„i and $m are eigenvalue and corresponding eigenvector of M^. Thus, we can calculate exact value of f/'C'^! k) 
by diagonalizing Mk and writing initial state of the photon in Fourier domain. The state of the photon after n steps 
in original domain "(/"(n, q) is given by the inverse Fourier transform of ^(n, fc), defined by Eq. (|A6|) . The eigenvalues 
of Mk are Ai = —1, A2 = 1, A3 = e"*"'' and A4 = e^^'' , where ujk £ [7r/4, 37r/4] and defined as coswfe = cos/c/v^. The 
corresponding eigenvectors are 



1 /3 + 2\/2cosfc 
' ''^2V 2 + V2cosfc 



V2+e- 
V2+e 



1 /3-2V2cosfc 
' ^' ~2V 2-\/2cosfc 



1 



1 



I $3) = 



2v/l + sin^fc 



l$4 



2v/l + sin2fc 



-l + \/2e-*('^'=--'') 

-1 

1 - x/2e"*('^^-+'') 

1 

-l + %/2e*("'=+'') ' 

-1 

1 - V2e'('^'="'') 

1 



If the initial state of the photon is |0, 0, h, x) the corresponding state in fourier domain will be ^"(0, k) = [1, 0, 0, 0]^ 
for all values of k. Similarly if the photon starts QRW from state |0,0, h,y) the initial state in Fourier domain will 
be [0,1,0,0]"^. Now from Ea. (|A9[) . we calculate the state of the quantum walker ip{n,k) in fourier domain. Taking 
the inverse Fourier transform of ip{n, k), we get the state of the walker after n steps il^{n, q) in the original coordinate 
space. For initial state of the photon |0, 0, h, x), we get 

\iphx{n,q)) = fhx\n,q,h,x) + fhy\n,q,h,y) + fy^\n,q,v,x) + fyy\n,q,v,y), 

where the coefficients fafs are given by 



(AlO) 



Jhx — 
Jhy — 
Jvx — 
Ivy = 



1 + (-1)"+" 


1 + (-l)"+9 


Stt 
1 + (-l)"+9 


87r 



dk 



^ COs{kq) 
7r 2 — \/2cos/i; 

"" V2cos[k{q+l)]-cos{kq) 
-TV 2 — V2cos/c 

"" C0s[fc(q + 2)] 



'^ 3 - 2V2cos(cjfe -k) _ 



1 + sin^ k 



i(nu)h+kq) ^^ 



'-dk 



dk 



-2ik 



Stt 



-TT 2 — \/2 COS k J_^ 1 + sin^ k 

"" V2cos[k{q + 1)] - cos[fc(g + 2)] 



-TT 1 + sin k 



+'"^Uk 



2 — \/2cosfc 



dk 



-k) 



1 



1 + sin^ k 



g~i(nc^,+fc9)^^ 



(All) 
(A12) 
(A13) 
(A14) 



Similarly, for initial state of the photon |0, 0, h, y), the state of the walker after n steps is 

li^hyin, q)) ^ghx\n,q,h,x) + ghy\n,q,h,y) + g^^\n,q,v,x) + g^y\n,q,v,y), 



where 



ghx 



9hy 



9vy 



1 + (-l)"+9 


Stt 




1 + (-!)"+« 


Stt 




1 + (-!)"+« 


Stt 




1 + (-!)' 


l+q 



" %/2cos[fc(g- l)]-cos(fcg) 
-TT 2 — \/2cos/i; 

'^ (3-2V2cos/c)cos(fcg) 



dfc 



"" 1 - V2e-*("'=-'=) _,-( 



1 + sin^ k 



nujk+kq)^j^ 



2 — \/2cosfc 



dfc 



TT -i{ni^k+kq) 



7r 1 + sin fc 



-dfc 



^ V2cos[k{q + 1)] - cos[fc(g + 2)] 



dk + 



^^-t{uk+k) _ I 



Stt 



-TT 2 — V2cosfc 7-7r 1 + sin^ fc 

'^ 2 cos fcg - 2V2 cos[A:(q + 1)] + cos[A:(g + 2)] 



g-»(nt^,+fcg)^^ 



2 — \/2cosA: 



^dfc 



IT g-i(n(^fc+feg) 



7r 1 + sin fc 



-dfc 
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(A15) 

(A16) 
(A17) 
(A18) 
(A19) 



Because of the factor 1 + ( — 1)"+'', coefficients fafs and ga/3 are nonzero for even values of n + g only, which is 
corresponding to the positions of detecting photon in the arrangement (see Fig[T]) for a fixed value of n. The integrals 
inside the bracket can not be evaluated exactly. Further the first integral inside the bracket is independent of n, and 
is responsible for the constant spikes in the probability of detecting photon near the initial position q = 0. The second 
integral inside the bracket completely characterizes QRW. 

In order to understand the nature of QRW we do asymptotic analysis of the second integral inside the bracket in 
the coefficients faf3 and ga[3 as follows. The integral has the form of 



/(n) 



f{k)e''"t'^'''^dk, 



(A20) 



where (/'(A;) = —j^k + fca), a — q/n, and f{k) is function of k only. In the limit of large n we use stationary phase 



-(wfe 

f25^ 



approximation [25^ and find the approximate value of /(n). 
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(A21) 



where fco ~ sin (— a/Vl — oi^) and prime denotes the derivative with respect to k. For a — ±l/\/2, <f){k) has 
stationary point of order two at A; = ±7r/2, i.e. 0'(±7r/2) = 0, 0"(±7r/2) = 0, (/)"'(±7r/2) ^ 0. At these points 



'(") " ('^)''' ^^^^/(-/2)e™^(--/^), for a = i=. 



I{n) 



ey^^ V2r(i/3) 

V3 



/(7r/2)e"'^(^/2\ for a 



1 

V2" 



(A22) 
(A23) 



Everywhere else /(n) has no stationary point and averages to zero for n — > cx). 
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